FLUCTUATIONS IN THE MIXED EVENT TECHNIQUE 
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ABSTRACT 

A method for an evaluation of fluctuations in the mixed event technique 



> 

£— n ' is proposed. It is shown, that, generally, the magnitude of the fluctuations 

00 ' 

£N| ' is proportional to A 3 ' 4 , where A is the number of produced events, which 

r->: 

^•-r ' should be compared with A 1 / 2 for the case of independent events. The 

On 

■"-^ ' formula eligible for the use in an analysis of experimental data is provided. 

q^ PACS number(s): 13.85.-t, 13.85.Hd, 25.70.Pq 



The mixed event technique was initially proposed [1] for the generation of 
uncorrected pion pairs for the intensity interferometry of identical bosons. 
Later on, the technique was widely used for the generation of background dis- 
tributions for different normalized correlation functions and in many other 
applications, e.g. for the calculation of background mass spectrum in reso- 
nance reconstruction. To use all the available statistic very often the same 



1 On leave from Moscow Engineering Physics Institute, Moscow, 115409, Russia 



particle is used in generation of many "events" ; thus produced "events" are 
not totally independent. Our goal is to evaluate the fluctuations in the pro- 
duced distributions. In this note we consider particularly the case of the 
distribution of pairs of particles in their relative momentum, but the method 
can be applied to any distribution generated by event mixing. 

Below we derive the formula for the variance of statistical fluctuation 
in the number of pairs in some region of the particles relative momentum. 
To generate pairs we use M pions from M different events (one pion from 
each event), generating altogether N%£ irs = M(M — l)/2 pairs, distributed 
among different bins in accordance with the value of relative momentum 
Qij = Pi — Pj] i,j = 1,2, ...,M. We study the fluctuations in one particular 
bin. It is convenient to introduce the function 

{1 if qij belongs to the bin 
(!) 
otherwise 

Using this notation the mean number of pairs in the bin under study 



N bin = N tot — _ prtot — (2) 

JV pairs JV pairs'^ — JV pairs"'- \*J 



Note that n depends only on the particle distribution in momentum and bin 
size and position, but does not depend on M\ 



N pairs i s a random variable (below we use the shorter notation X = Ni^ rs ) 
with a mean value X and variance ax; 

To get the variance we have to calculate 



a\ = X^-X 2 



- 4l>,) 2 -(^r^)V = 

= i{M(M-l)(M-2)(M-3)(n-n^-n 2 ) + 
+4M(M - 1)(M - 2){n-^i - n 2 ) + 
+2M(M-l)K-n 2 )}. (4) 

In the first term of the last expression all four indices (i,j,k,l) are dif- 
ferent; consequently njjrTki = n 2 . Due to this the first term equals to zero 
and does not contribute to the variance. Let us go now to the last term. Its 
contribution is 

^1^%-^xa-w), ( 5) 



where we used the equality n 2 k = n, following from the definition ([]]). Usually 
n<l, thus the contribution of the last term to the variance is just the mean 
number of pairs in the bin; the value which we would expect if the pairs 



were totally independent. But it is known that the pairs are really correlated 
because each pion is used in generation of (M — 1) pairs. The price for 
this is the second term in the expression for the variance. Let us denote by 
a = no-Tiki — n 2 . Note that this value, like n, does not depend on M. Then 
for the variance we have: 

o\ = M(M - 1)(M - 2)a - X(l - n^) « M 3 a - X. (6) 

If aM 3> n the first term becomes dominant and 

a\ w a(2/nf 2 X 3/2 (7) 

The fluctuations in this case are proportional to (^pa" rs ) 3//4 , and can be much 
larger than for the case of totally independent pairsQ. 

The parameters a and n can be expressed using one-particle density: 

1 f 
n= =2 / dpidpjp^ip^p^ip^riif, (8) 

N J 
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a = n 2 + =3 / dpidpidp k p {l \pi)p {1) ipi)p^{p k )n ik n k i, (9) 



where N is the mean number of particles (the value one-particle density is 

normalized to). We used the fact that the individual particles in our tech- 
2 This is exactly the law seen by Zajc et al. [2] empirically. In the Ref.[2] this law was 

derived in a simple model, but for the real case the authors studied empirical relations. 



nique are independent, and we can express two- and three-particle densities 
as a product of one-particle densities. 

In the analysis of the experimental data these parameters should be sub- 
stituted by their estimators (denoted by ()) directly using the data. For 
example, one can use 

ns! (riii) = rrTTT ^y^nij (10) 

x 3l M(M - 1) j£ J y ' 

a « (n ik n kl ) - (n^) 2 = M , M _ 1 \, M _ 2 ) ^ nikUkl ~ ^'^' ( U ) 

where in the last sum all three indices (i, k, I) are different. Since a does not 
depend on M, the average does not need to be carried out over the entire 
sample. 

Once we have seen that the fluctuations in the number of mixed pairs can 
be rather large it becomes worthwhile to recall that the estimator usually used 
for the correlation function defined as a ratio of estimators for the numbers 
of the true and mixed pairs is biased: 



T\jtrue ivrtrue 

-T^T = ^pavrs , c true = jVilL ( i2 ) 

\ z / jymixed I 2, mixed v ; 

pairs ^-^ pairs 



The necessary correction can be derived from the equality 



1 1 I., (X-X) 2 (X-X) 3 . 1 .„ (ri. 

x x + (x-x) x { + x 2 X 3 j x l x 2J ' 

(13) 
where a\ is defined by Eq.6. The corrected estimator for the correlation 
function is 

2 

(C 2 ) corrected =(C 2 )(l-^), (14) 

where, as earlier, we denote by X the mean value of mixed pairs in the bin 
under study 
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